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Abstract 

In this work we present a solution of the Boussinesq equation. The derived formulas include solitons, 
Schwartz class solutions and solutions, possessing singularities on a closed set Z of R 2 {(x,t) domain), 
obtained from the zeros of the tau function. The idea for solving the Boussinesq equation is identi- 
cal to the (unified) idea of solving the KdV and the evolutionary NLS equations: we use a theory of 
evolutionary vessels. But a more powerful theory of non-symmetric evolutionary vessels is presented, 
inserting flexibility into the construction and allowing to deal with complex- valued solutions. A powerful 
scattering theory of Deift-Tomei-Trubowitz for a three dimensional operator, which is used to solve the 
Boussinesq equation, fits into our setting only in a particular case. On the other hand, we create a much 
wider class of solutions of the Boussinesq equation with singularities on a closed set Z. 
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1 Introduction 



The Boussinesq equation [Bou72pl 



i:) 2 



qtt = nsl&I** ~ W] W 



is a foundation of a shallow water theory. A fundamental solution of this equation was presented by P. 
Deift, C. Tomei and E. Trubowitz in [DTT82] , The basic idea in this work, presented by V. E. Zakharov 
[Zac74| . is a development of a scattering theory of a three dimensional operator 

~ _ . d 3 I d d 
dx 3 i dx dx 

d 2 

and using the fact that for Q = i(3 — — — 4q), the operators L, Q constitute a Lax pair: 

dx 1 

^-L = [Q, L] = QL - LQ. 
dt 

It is assumed in |DTT82) that p(x),q(x) are in a Schwartz space, but the theory goes through for 
q(x),p(x) with only a finite number of derivatives and a finite order of decay. The scattering data in this 
case consists of a list of 6 functions, whose evolving with t under an analogue of Q is studied producing a 
solution of the Boussinesq equation |T} with a given initial value q(x, 0) = q(x). Since we will use instead 
of L its multiplication on — i, we define 

L = - ll = i- 2q l~^ +ip) (2) 

and a very particular (inverse) scattering of this operator will be researched. More precisely, we will 
discuss solutions of0 



In this setting q(x) is an "arbitrary" function, but p(x) is derived from q(x) (see formula (|38[) ) using the 
formula p'(x) = P(f q, q, q 1 , . . . , q^ n ') for a polynomial P, which is not derived explicitly, since we do not 
use it. 

In order to solve Q, we present a similar to [DTT82] scheme, where we use a special case of the 
inverse scattering theory: we create inverse scattering of p]l. where p(x) is not arbitrary and then evolve 
q(x) with t. Let us explain first the (inverse) scattering theory. The scattering data is encoded in our 
setting in a matrix-valued function. In fact, there is an almost "one-to-one" correspondence (Theorem 
I14[) between the coefficients q(x), defining L ((2| and 3x3 matrix- functions S(X) possessing a realization 
[BGR90] 



S(X) = I- C Xo '(A/ - Ay 1 B 



1 

1 

1 



AX + X A C + BoarCo = 0. 



Q2 

1 The full Boussinesq equation qtt = — ~[q + qxx — 4q 2 ] was shown by McKean lMcK78| to be equivalent to (flj, significantly 

dx 2 

simplifying algebra computations. 

2 we will usually denote by q' the partial derivative — q, and by q the partial derivative — q. Similarly, for the higher 

dx " dt 

d 2 

derivatives: q" stands for — -q, etc. 

dx 2 
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Here for an auxiliary Hilbert space H the linear operators act as follows: Co : T-L — ¥ C , Ac; ,Xo, A : T-l — ► T-L, 
Bo : C 3 — > rl. Let us assume for the simplicity of the presentation, that all the operators are bounded. In 

, 1 

order to construct q(x), uniquely defined from S(X), use the following 4 steps, fixing 01 = 











-1 



solve for B(x) from B'ai = —ABo 2 — By (J3J) with initial B(xo) = Bo, 
solve for C{x) from o\C' = —oiCA^ + Cy ((6} with initial C(xo) = Co, 
solve for X(x) from X'(x) = B(x)a 2 C(x) Q, X(x ) = X , 

define y*(x) — 7 + a 2 C (x)X~ x (x) B (x)ai — aiC(x)X~ 1 (x)B(x)a 2 ©, for all points where X(x) is 
invertible. 

Then we prove in Theorem [4] that the function 

- C(x)X~ 1 (x)(XI - Ay 1 B(x)ai 



S(X,x) 



d 3 



is a Backlund transformation for the operator L from the trivial Lq = — - to a more complicated one 

dx s 

3 d 2 

L ([2]), in which q(x) = ~2~d~~2 m< iet (^o~ ^-( x )) an< i p( x ) i s defined from q(x) up to a constant. In the 
regular case, explained here the coefficients q(x),p(x) are analytic functions at all points, where X(x) is 
invertible. 

Letting the operators C, B, X further evolve with respect to t, we will obtain a solution of the Boussi- 
nesq equation ([T}. In order to show it, we take three additional matrices 



(71 = (71, 



(72 = 



" 


— i 


" 




" 





" 


I 








» 7 = 











_ 
















i 



B{x), 



and follow these steps (last step is the same as the previous fourth step): 

4. solve for B(x,t) from Bo\ — —ABa 2 — By (|5Up (see footnote with initial B(x,to) 

5. solve for C(x,t) from o\C = —a-iCA^ + Cy (151) with initial C(x,to) = C(x), 

6. solve for X(x,t) from X(x) = B(x,t)a 2 C'{x,t) ©, X(x,t ) = X(x), 

7-4. define 7,(2, t) = 7 + a 2 C(x,t)X~ 1 (x,t)B(x,t)a 1 - oiC{x, t)X' 1 (x, t)B(x, t)a 2 ©, for all points 
where X(x, i) is invertible. 

3 9 2 



It turns out that q(x, t) = — — — j lndet ( 



( ; , y X X(a;,i)) (f3T)l and satisfies the Boussinesq equation fT} (see 

Theorem 1 1 71 for details). Construction of the coefficient q(x) from a realized function S(X) is called </ie 
standard construction of a vessel and is presented in Section 

The simplicity and richness of this construction is best revealed in soliton formulas (Section I6,I[1 . 
By choosing the inner space H = C - the one dimensional Hilbert space, we create a classical soliton 
9fi 2 



q(x,t) = 



2,V3 



58} and another one q(x, t) = 



2 cosh (---/i(a; + f/^)) 







3 2\/3:e/j 



+ e 



59)1 . Here /i is 



an arbitrary complex parameter. 

This construction seems to be more concrete and suitable for the solution of the Boussinesq equa- 
tion fTJ), because it uses just "enough" of the very powerful and complicated inverse scattering theory, 
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developed in |DTT82j . Notice that in this later work the coefficients q(x),p(x) are quite arbitrary. In 
our work, on the other hand, the coefficient p(x) is uniquely determined (up to a constant) from q(x). 
Still, the formulas enable to produce solutions, applying the standard construction of a vessel, explained 
above, to different S(X). Moreover, at the same section [3] we show that much more general classes of 
solutions arise, as we impose as few as possible restrictions on S(X). The solutions, presented in [DTT82] 
are either in the Schwartz class or exponentially decaying and correspond in our setting to an analytic 
5(A), possessing jumps along the real negative axis. 

Finally, this work presents a very general setting for a construction of solutions of ((TJ). We find 
necessary regularity assumptions (|61[) on the operator B(x, i) such that we can create a vessel, and hence 
a solution of {TJ. This is proved in Theorem 1 191 Following the remark after Theorem 1 191 one can easily 
construction a solution, which fails to be five times a;-differentiable (four times are necessary for the 
existence of |T|)). 

2 Scattering theory of the operator L 

We start from the definition of the vessel parameters, which create an inverse scattering theory of L ([2]). 
Definition 1. The vessel parameters are defined as follows 





' 





1 " 




1 





~ 




" 








0~1 = 





1 





,0-2 = 











,7 = 








1 




1 
























-1 






2.1 Non-symmetric regular vessel, realizing a scattering theory of L 

Definition 2. A (regular, non-symmetric) vessel, associated to vessel parameters (see Definition^ is a 
collection of operators, spaces and an interval I 

9Jr es = (C(x),A c ,X(x),A,B(x)-a 1 ,a 2 , 1 ,- ft (x)-H,C 3 ;V), (4) 

where the bounded operators C(x) : H — > C 3 , A(,~K(x),A : H — > H, B(x) : C 3 — ¥ H and a 3 x 3 matrix 
function 7* (a;) satisfy the following vessel conditions: 

= -(AB(T2 + B^ 1 , (5) 
= <T^ 1 ( 7 C-a 2 CA ( ), (6) 

= Ba 2 C, (7) 

7* = 7 + aaCX^Bo-i -a 1 CX~ 1 Ba 2 , (8) 

iX + Xic =-B(7iC. (9) 

The operator X(x) is assumed to be invertible on the interval I. If Aq = A* and C = B* we call such a 
vessel symmetric. 

Remarks: 1. Notice that the operators C(x),X(x), B(x) are globally defined for all C(x), B(x) 

are solutions of operator-valued differential equations with constant coefficients and X(a;) is obtained from 
them by a simple integration. 2. For the definition of the matrix-function 7* (2) we need the invertability 
of the operator X(x), so in general we may suppose that 7* (a;) is defined for all x £ R, except for those 
points where X(x) is not invertible. For simplicity, we take an interval I and in a more general setting 
we consider 7,(2;) on R, except fro the points where X(i) is not invertible. 



dx^ 
dx^ 
dx 
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Theorem 3 (permanency of the Lyapunov equation). Suppose that B(x),C(x),X(x) satisfy ([5]), ((6} and 
(J7J) respectively, then if the Lyapunov equation @ 

AX(x) + X(x)A ( + B{x)a 1 C(x) = 

holds for a fixed xq £ I, then it holds for all x el. 

Proof: By differentiating the left hand side of ([9]), we will obtain that it is zero. □ 
By the definition, the transfer function of the vessel QJ re9 is defined as follows: 

S(X, x)=I- C(x)X~ 1 (x)(XI - Ay 1 B(x)a 1 . (10) 

Notice that poles and singularities of S with respect to A are determined by A only. We would like to show 
that the function S(X, x) realizes a Backlund transformation of the corresponding LDEs: multiplication 
by the function S(X, x) maps LsOl, AMV12, Melll a solution of the input Linear Differential Equation 
(LDE) with the spectral parameter A 

A<72it(A, x) — oi-^-u(X, x) + 7«(A, x) = 0. (11) 

to a solution of the output LDE with the same spectral parameter 

d 

Xa 2 y{X,x) - ex— y(X,x) + j*(x)y(X,x) = 0. (12) 

The function 7, is defined by the Linkage condition ((8|. The fundamental solutions of (|11[1 and (1121) . 
which are equal to / (the identity matrix matrix at x — 0) are denoted usually by $(A,:r) and $»(A,x). 
In other words, these are matrix functions satisfying: 

\o 2 HKx)~a 1 -?-$(\,x)+ 1 $(\,x) = 0, $(A,0)=J, (13) 

\o- 2 $,(\,x)-o- 1 -^$,(\,x)+~/,(x)$4\,x) = 0, *»(A,0)=7. (14) 

Theorem 4 (Vessel=Ba,cklund transformation LsOl AMV12, Melll J. Suppose that u(X,x) satisfies 
(|11[) . then y(X,x) = S(X, x)u(X, x) satisfies (1121) . 

Remark: Notice that it is enough to prove that the transfer function satisfies the following ODE: 

^S(A, x) = o^(a 2 X + f*(x))S{X, x) - S{X, x)o^(o 2 X + 7). (15) 

Proof: First we calculate 

^ (C(x)X- 1 (x)) = o--\jC - OiCA^X' 1 - CX^B^CX' 1 = by © 

= a-^CX- 1 + o- 1 a 2 CX- 1 A + o-- 1 o- 2 CX- 1 Bo- 1 CX~ 1 - CX- 1 Ba 2 CX~ 1 
= ai 1 a 2 CX~ 1 A + ^(7 + o 2 CX- 1 Ba 1 - aiCX" 1 Ba^CX' 1 = by © 
= o- ] : 1 a 2 C*X- 1 J 4 + f 7- 1 7»CX- 1 . 



So we obtain that 



d (CX- 1 ) = al x a 2 CX~ x A + arS-CX -1 . (16) 



dx 
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Let us differentiate next the transfer function using |10 



4-S(X,x) = —^(CX- 1 )(XI-A)- 1 Bai-CX- 1 (XI-A)- 1 -^Bai = by ©, CE1 
dx dec dx 

= (a- 1 a 2 CX- 1 A + a-^/^CX- 1 ) (XI - A)~ 1 Ba 1 - CX _1 (A/ - A)~ l (A Bo 2 + B-y) 

= a^^ls - J) - (5- iVrS + a 1 - 1 cr 2 CX- 1 A(Al - A) _1 B(7i - CX _1 (A1 - AJ^ABua 
= insert A = A ± XI and expand 

= a^y^S-I) - (S - I)(T^~t + a- 1 a 2 CX- 1 Box - Xa^ 1 cr^CXr 1 {XI - A^Bai- 

+CX~ 1 Ba 2 - ACX -1 (A7 - A)- 1 Bo 2 
= using ©, JTDJ 

= ^.(s - /) - (5 - /yrS + o-i _1 (7* - 7) + Ao-rVaS - so-r^ 

= (T 1 " 1 (Aa2 + 7»)5 — Serf (A<X2 + 7). 

□ 

Expanding the transfer function S(X,x) into a Taylour series around A = 00, we obtain a notion of 
the moment: 

H n (x)ai 



S(X,x)=I-^2 



where by the definition the n-th moment H n (x) of the function S(X,x) is 

H n (x) = C(x)X" 1 (x)A n B{x). (17) 

Using the zero moment, for example, we obtain that the so called "linkage condition" ([8]) is equivalent to 

7*(x) = 7 + a 2 H (x)a! - o- 1 H {x)a 2 . 

There is also a recurrent relation between the moments H n (x), arising from (|15l) : 
Theorem 5. The following recurrent relation between the moments of the vessel QJ r e 9 holds 

a^ 1 a 2 H n+1 - Hn+io-'^i 1 = (H n )' x — cr^ + Hnja^ 1 . (18) 

H n (x)oi 



Proof: Follows from the differential equation (|15l) by plugging S(X, x) — I — ~}2^L 
Let us investigate more carefully the LDEs (|lll) and (|12[1 . Denote u 



Wl 
«3 



\n+l 

then (|11|) becomes 



□ 



Ami 




^3 









" " 







u' 2 


+ 














1 




-U2 








Solving this we obtain that 

( u 2 — —u[, 

< u 3 — u 2 = (19) 
^ u'i = —Aiii. 

We can see that actually this equation is equivalent to a third-order differential equation with the spectral 
parameter A: 

u"' = -Xui. (20) 



G 



In order to analyze (|12p . we denote first moment Ho(x) = [tt^] 
the linkage condition ((8| becomes 



TTll 


7T12 


7T13 


7T21 


7T22 


7T23 


7T31 


7T32 


7T33 



, and as a result, 



7* =7 + 



7T13 — 7T31 7T12 Till 
-7T21 

-7T11 



Denote next y = 



V2 
2/3 



and plugging the expression for 7* just derived into ()12[). we will obtain that 



Xyi 




2/3 




(nvs — n-si)yi + 7T12J/2 + TTlll/3 




" " 







2/2 


+ 


— 7i"2iyi + 2/3 
















-TTnyi — y 2 








or solving this: 



J/2 = -7TUJ/1 - £/!, 

J/3 = 7T 2 iyi + J/2 = 7T21J/1 - J/l^"ll - KllVl 

y'i - Zqy'i - (<?' +p)yi = -Xyi, 



2/1: 



where 



I'll + 7r 12 + ""21 — 27TX! 



p(x) 



1"12 — 7T21 • 



— i(— 7!"13 + 7T31 + 7Tll(7Tl2 — 7T21 ) 

In other words the equation (|12[1 is equivalent to 

2/1 ' - 2gyi - (g' + ip)yi = -Xyi. 

Now we are ready to justify the term "scattering matrix" attached to S(X, 0). An independent set 

" 2/i (A, x) 
y(X,x)= yi(x,x) 
V3(X,x) 

of solutions of can be derived from (|12|l in the following form 



(21) 

(22) 
(23) 





' 1 " 




y(X,x) = S(X, x)$(X, x) 





( 












1 " 




$,(X,x)S(X,0) 





) 










The fundamental matrix "1?(A, x), solving (|lip with the initial condition "^(A, 0) = / (-the identity matrix) 
is obtained from (|19[1 and can be explicitly written as 

--R3 - 



i-2 

^2 



fc 
fll 



fc7?3 #1 
-fc 2 i?2 -fc/Js 
7?1 =Q 2 (Si+ J B 2 +S3), 

i? 2 =a 2 £i + a£ 2 +-B3, 
i? 3 = a 2 Ei +E 2 + aE z 



(24) 



where E\ — e 



,E 2 



,Ea 



for a 



,27ri/3 / 3 



(a = 1). This matrix is analytic in A, because 



examining Taylor series of Ei we will come to the conclusion that all the entries of $ depend on fc = A. 
The structure of S(X, x) is also known from (|10[) . so we can study solutions of ([2| or equivalently of (|12|) . 
creating in this manner the (inverse) scattering of L @. 
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2.2 Structure of the moment Hq(x) 



Let us examine the recurrence relation (|18[) . We will research for the simplicity of the presentation the 
structure of the first moment Ho(x), but almost the same structure will actually apply for all moments. 
Let us denote 



HUx) = 



9u 


912 


913 


921 




923 


931 


932 


933 



Then the left hand side of (1181) becomes 



Ci 1 <J2Hi(x) — H\(x)o2<J x 1 = 



-gn 
-gal 
3ii 312 3i3 - 33i 



which must be equal to 



(H )' x - a^H + Hoja- 1 = 

Til + ""12 + 7T21 + ir'n 
7T117T21 — 7T31 + 7T22 + 7T 2 i 
— 7T117T13 — 7T127T21 + 7T32 + 7T 31 



7T117T12 — 7I"13 + 71~22 + n' 12 
7T127T21 — 7T23 — 7T32 + 71~22 

Res(32) 



7T117T13 + 7T23 + Tt'l3 

Res{23) 
Res(33) 



where 



7?es(23) = 7T137T21 — 7T33 + 7r 23 , 

Res(32) = — 7Tl2(7Tl3 — 7T31 + 7r 22 ) — 7T117T32 — 7T33 + 7132, 
7?es(33) = — 7T127T23 + VTl 3 ( — 7Tl3 + 7T 3 l) — 7111^33 + 7I33 . 



Equating we obtain that 

311 = — 7T117T13 — 7I"23 — ^13i 921 

and the following system of equations: 



-Res(23), 312 = -Res(32), 313 — 331 = i?es(33) 



(25) 



TTll + 7T12 + 7T21 + 7Ij.i = 0, 
7T117T12 — 7T13 + 7T22 + 7r 12 = 0, 
7T117T21 — 7T 3 1 + 7T22 + 7T 21 = 0, 
7Tl27I"21 — (tT23 + 7T32) + 7I~22 = 0, 
— 7TH7T13 — 7T23 — 7T 13 = — 7Tll7ri3 



- 7T127T21 + 7T32 + 7T 31 (= 311) 



or rearranging 



(7T12 + 7T1! 



7T21 

7T22 = — 7Tll7ri2 
7T22 = — 7T117T21 



TllJ, 
7T13 — 7i~12, 
""31 - 7T 2 1, 



(26) 



7T32 = 7T12 7T21 — 7T23 + ^22 j 
7T32 = 7T127T21 — 7T23 — 7T 31 — Tl' ls . 

Plugging the fourth equation of (|26|) into the last one, we will obtain that the later becomes n' 31 + 7r 13 = 
— 7r 2 2 or requiring a normalization 



we obtain 



7T3i(a;o) +7ri3(a;o) = ~K22{xt>) 



7T31 + 7T13 = — 7T22- 



(27) 
(28) 
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In a similar manner, one can solve for some other entires and we summarize these intermediate results 
in the next lemma. Additional relations between the entries of Ho(x) are obtained if we consider (|18[) 
for n = 1. For example, formulas similar to (|26p are as follows: 

Tiiflii + 312 + 92i + g'n = 0, 

7T11312 - 313 + 322 + 9'l2 = 0, 
ffLl7T21 - 331 + 322 + 321 = 0, 
312 7121 ~ (323 + 332) + 322 = 0, 

— Kligi3 - 323 - 3l3 = — THS31 + 3ll(7r31 - 7Tl 3 ) - 71i 2 #21 + 332 + 331- 

From these formulas follow the folllowing relations 

+ flia + 321 + 3u = 0, , 2g . 

7T11312 - 7Tllfl 2 l - (ff!3 - 33l) +ff'l2 -321 = 0, 

which can be rewritten as relations on the entries of Ho(x) in view of (|25[) . Other three relations, which 
are obtained are 

f 7T11312 - 313 + 322 + 312 + 3ll1"21 - 331 + 322 + 321 = 0, 
< 312 7121 ~ (323 + 332) + 322 = 0, 

[ -7T11313 - 323 - g'li = -7T11331 + 311(7131 - 7Tl 3 ) - 71l 2 <721 + 332 + 331, 

which serve to find g 22 , 323+332, 313+331- It turns out that actually the relations (|25|) are not independent 
and a formula for 7112711! is derived from them. All these results are summarized in the next Lemma and 
we notice that the exact formulas for 7i 23 , 7i 33 are omitted, because we are not interested in their form: 

Lemma 6. The following relations between the entries of the first moment Ho(x) hold 

7121 = — (7112 + 7tii + ir'n), (30) 

7122 = — 71117112 + 7T13 — 7lJ 2 , (31) 

7131 = 7T13 — TTi! — 7Tn(27ri2 + 37rn) — 27r 12 — 7T1!, (32) 

7132 = 71127121 — 7123 — 711271^ — -((ttJi) 2 — 71llTTii) + -71^", (33) 

lis = -^Wi) 2 + 7Ui7ri 2 + ^TiiiTi^ + < 2 + iyri'l, (34) 

6711271^ =-(67111 + 1571^)^1+3711171^+71^. (35) 

The relations (|29f) result in formulas for 7i 23 , 7i 33 . 

Proof: Notice that we have obtained 5 relations on the entries of Ho(x) in (|26|) . 2 relations in (|29p and 



one more relation from the fact that (I25|) together with (|29p are overdetermined. It is possible to derive 
these formulas using any symbolic computation software. □ 
It would be interesting to derive formulas for the entries of Ho(x), relying on the structure of the 
matrices <7i,<72,7 only, without referring to tedious direct computations. It can be probably done using 
the following notion, appearing first in M elll] : 

Definition 7. Tau-function of the vessel QJ is defined as follows: 

t(x) = det(X _1 (xo)X(x)), (36) 

where xo £ I is an arbitrary point. 

The fact that this object is well-defined follows from an equivalent to (0 equation 



X(x) =X + B(y)a 2 C(y)dy, 
Jo 

so that X ( ^ 1 X(x) — I + T(x) for a trace-class operator T. Using this notion we obtain the following 
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Theorem 8. The coefficient q(x) possesses the following formula: 

q{x) =-l£ lnT ^ = -|^ 7ri1 - (37) 
The derivative of the coefficient p(x) is a differential polynomial in tyu: 

p (x) = P(wn, tt^i, . . .), P - a polynomial. (38) 

Proof: Inserting (|30|l into the formula (|22|l for q(x) we find that 

, s 7Tll + 7T12 + 7T21 — 27TX! 3 , 

= ^ = - 5 7T U . 

Using a formula for the determinant of an operator, we obtain that (sec GK69, AMV12 for details) 
t'(t) 

— = tr(X'(x)X _1 (2;)) = tr(B(x)a 2 B*(a;)X" 1 (a;)) = tr(a 2 H (x)) = n- L1 
t(x) 

and the result follows for q(x). As for p(x), we differentiate the formula for p(x), appearing in (|22l) : 

p (x) = 7Ti 3 +7T 3 i) - « — (7Tn(7ri2 - 7T 2 l) - )• 

Then using (|32|l . (|34|) . (|30[) and (|35f) we will obtain a differential polynomial in 7Tn. □ 
One can also derive the following formulas, corresponding to the symmetric case C — B* , Aq = A* , 
plugging the definition (|36[1 and using the fact that Ho(x) is self-adjoint: 

Lemma 9. For the symmetric case, C = B* , A( = A* , the following relations between the entries of the 
first moment Ho(x) = Hq(x) hold 

It It It 

5K7Tl2=-- , 7T22 = , 5R7T13 = , 

It 3 t 6 t 

£S7T13 = 7Tll$j7ri2 + 5*71x2, 
^TT 2 9 1 1 T< 4 > 

K7r2 3 = ^ + 2 (g(x) 2_l ?>))+ Xr_ 

3 Vessels with unbounded operators. Standard construc- 
tion of a vessel 

The ideas presented in this Section can be found in [Mela] for the symmetric case. The class of functions 
serving as "initial conditions" for the transfer functions of vessels is defined as follows 

Definition 10. Class 7t(o~i) consist of p x p matrix-valued functions 5(A) of the complex variable X, 
possessing the following representation: 

5(A) = I - CoX" \XI - A)~ 1 B o 1 (39) 

where for an auxiliary Hilbert space H there are defined operators Co : "H — > C 3 , Ac, , Xo, A : H — > H, 
Bo : C 3 — > H. A general matrix-function 5(A), representable in such a form is called realized. Moreover, 
the operators are subject to the following assumptions: 

1. the operators A,A^ have dense domains D(A), D(A() . A, A^ are generators of Co semi-groups on 
H. Denote the resolvents as follows R(X) = (XI — A)^ 1 , R((X) = (XI + Aq)^ 1 , 

2. the operator Bo satisfies R(X)Boe 6 H for all X £ spec(A), e 6 C 3 , 
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3. the operator Xo is bounded and invertible, 

4- the Lyapunov equation holds for all A £ spec( J 4) U spec(— A( ): 

R(X)Ko - XoJic(A) + R(X)B (tiCoR C (X) = 0. (40) 

We call an element of 7Z.((Ti) as scattering matrix- function. The subclass li(ai) C 7?.(<ti) consists of 
symmetric functions, i.e. satisfying S(X)a^[ 1 S* (— A*) = o-j" 1 . The Schur class SU(ai) C U consists 
of symmetric functions, for which Xo is a positive operator. The sub-classes of rational functions in 
SU,U,'R. are denoted by rSU,rtl,r'R, respectively. 

When 5(A) is just analytic at infinity (hence A must be bounded), there is a very well known 
theory of realizations developed in [BGR90] . For analytic at infinity and symmetric, i.e. satisfying 
S* (— \)o-\S(\) = o\, functions there exists a good realization theory using Krein spaces (H is a Krein 
space), developed in JDLdSpl . Such a realization is then translated into a function in U(o\). The sub- 
classes U, SU appear a lot in the literature and correspond to the symmetric case. We will not consider 
these two classes in this work and refer to A MV12] , 

Now we present the standard construction of a vessel 03 from the given realized matrix-function 5(A) 

1. Let B(x) be the unique solution of the following equation (-R(A) = (AJ — ^4) _1 ) 

^i?(A)B = -(AR(X)Ba 2 + R(X)B-y)a^ 1 (41) 

satisfying B(xo) = Bq. This equation is solvable because the coefficients 0-1,02,7 are constant 
and A is a generator of a Co semigroup. For this equation to hold we must require the following 
regularity assumptions: 

VA £ spec(A) : R{X)B(x)o 2 C i C D(A), (42) 
VA £ spec(A) : R(X)B(x)>yC 3 C H, (43) 

2. Let C(x) be the unique solution of (|6]) with the initial condition C(xo) = Co, when we consider the 
equation ©, applied to vectors in the dense set D(A^) only, 

3. Solve (0 for X(x), satisfying X(xo) = Xo and let I be an interval, including xq on which %.{x) is 
invertible. Notice that B(x)a 2 C 3 C H follows from J42J. 

4. Define 7*(x) on I by (J8$. 

The main reason, why we call 5(A) as the "scattering data" is the fact that 7*(x) (generalized potential) 
is uniquely determined from 5(A) by this construction. The question of uniqueness of 5(A) for a given 
potential ~y*(x) will be studied in Section U A vessel in a more general form is defined as follows 

Definition 11. The collection of operators and spaces 

21 = (C(x) , A c ,X(x),A, B{x) ■ a! , era, 7, 7* (x) ; H, C p ; I), (44) 

is called a vessel, if C(x) : C 3 — > H, A$, X(x), A : H — > H, B(x) : C p — > H are differ entiable linear 
operators, subject to the following conditions: 

1. the operators A, Aq have dense domains D(A), D(Aq) . A, Aq are generators of Co semi-groups on 
H. Denote the resolvents as follows R(X) = (XI — A)^ 1 , Rc,{X) — (XI + A^) -1 , 

2. B(x) satisfies regularity assumptions (|42f> . (|43[1 . and the equation (|41|l . 

3 At the paper |DLdSI a similar result is proved for functions symmetric with respect to the unit circle, but it can be translated 
using Calley transform into S*(— X)a\S(X) = <J\ and was done in _Mclb, AM V12] 
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3. C(x) satisfies © on D{A^), 

4- X(x) is bounded, and invertible on I and satisfies 0, 

5. the Lyapunov equation holds for all x £ I, A $ spec(yl) U spec(— Aq): 

R(X)X(x) - X(x)R c (X) + R(X)B(x)a 1 C(x)R < (X) = 0. (45) 

6. 7* (a;) satisfies ([8]). 

The class of the transfer functions is defined as follows 

Definition 12. CZass I = %(o~i, o~\, o~2, 7; I) consist of 3 x 3 matrix-valued (transfer) functions S(X,x) 
of the complex variable X and i£l /or an interval I = [a, 6] possessing the following representation: 

S(X,x) = I — C(x)'%r 1 (x)(XI — A)" 1 B(x)ai (46) 

where for an auxiliary Hilbert space H, the operators C(x) : H —¥ C 3 , A^,X(x), A : H — > H and 
B(x) : C — > H constitute a vessel 27 (|44|) for some 0-2,7. 

An analogues of the Backlund transformation Theorem [3] for this new setting can be found in Mela 
for the symmetric case. Actually, the assumptions on the operators were found in such a manner that 
this theorem still holds. We omit its proof, since we are interested in solutions of {TJ. 



4 "Uniqueness" of the scattering data 



Let us consider now the uniqueness of the scattering matrix S'^O). First we prove the following 
Lemma 13. Suppose that we are given a regular vessel (f3| 

9Jre S = (C(x),A,X(x),A < ,B(x)-a 1 ,a 2 , 1 ,'y t (x)-H,C 3 -J), 

realizing coefficients q(x),p(x). Let S(X,x) be its transfer function, defined in (|10|) . Let Y(X) be an 
arbitrary 3x3 matrix function, commuting with the fundamental solution <!>( A, x) of 1)1111 . Then S(X, x) = 
S(X,x)Y(X) is the transfer function of a vessel realizing the same coefficients q(x),p(x). 

Proof: By the definition it follows that 



S(X,x) = $ t (X,x)S(X,0)$~ 1 (X,x). 



So, 



S(X,x) = S(X,x)Y(X) = $4X,x)S(X,0)$~ 1 {X,x)Y(X) = $»(A,x)S , (A,0)y(A)$" 1 (A,x) 

and realizes the same coefficients q(x),p(x). By the standard construction, there is a vessel 93, whose 
transfer functions is S^A, x). □ 
Let us investigate the structure of a matrix Y(X), commuting with $(A,a;). Using the form (|24|l . it is 
easy to conclude that a matrix, which commutes with $(A, x) must be of the form 



F(A) = a(A)/ + 6(A) 



1 


-A 











+ c(A) 



1 
-A 
A 



(47) 



by considering the coefficients of pure exponents in <&(X, x)Y(X) = Y(X)&(X, x). 

Theorem 14 (Uniquness of the scattering matrix). Suppose that S{X,x), S(X,x) are the transfer func- 
tions of two regular vessels 93 reg ,93 re9 , defined in Q, realizing the same potential 7*(x). Then there 
exists a matrix Y(X) G TZ. such that 

S(X,x) = S(X,x)Y(X). 
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Proof: Let us consider the function S 1 (A, x)S(X, x). By the definition this functions maps solutions of 
the input LDE {11} to itself: 

S _1 (A,»)S(A,x) = ($ t (\,x)S(\,0)$~ 1 (\,x)y 1 $*(\,x)S(\,0)$-\\,x) = 

Plug here, the formula (|24p and find conditions so that the coefficients of the exponents e~ kx , e~ akx , e~ a kx 
vanish. This is necessary for making this function bounded at infinity out of the spectrum of A. Then 
calculations show that actually S^ 1 (A, 0)S f (A, 0) must commute with $(A, x) so that this functions cancels 
all the singularities at infinity. As a result, by the preceding arguments it must be a function Y(A) of 
the form ((47)) . And we obtain that 

S- 1 (A,0)S(A,0)=y(A), 

from where the result follows. □ 
Another, weaker form of the uniqueness is used later in the text and is presented in the next Lemma. 
We emphasize that a similar theorem lemma was proved in the Sturm-Liouville case in [Melb] and in 
|Fad74] for purely continuous spectrum. 

Lemma 15. Suppose that two functions S(X,x), S(X,x) are in class X(<7i, C2,7), possessing the same 
initial value 

S(X,0) = S(A,0) 

and are bounded at a neighborhood of infinity, with a limit value I there. Then the corresponding outer 
potentials are equal: 

7*00 = 7*00- 

Proof: Suppose that 

S(X, x) = <f>*(X,x)S(X,0)$~ 1 (X,x), S(X,x) = $*(X,x)S(X,0)$- 1 (X,x). 

Then 

S~ 1 (X,x)S(X,x) = $,(\,x)$~ l (\,x) 

is entire (the singularities appear in S(X, 0) = S(X, 0) only and are cancelled) and equal to / (- the identity 
matrix) at infinity. By a Liouville theorem, it is a constant function, namely I. So $»(A, x)^~ 1 (X, x) = I 
or 

$.(A,aO = $,(A,aO. 

If we differentiate this, we obtain that 7*00 = 7*0)- ^ 



5 Choice of the parameters realizing the Boussinesq equa- 
tion on) 

5.1 Realizing the Boussinesq equation by a Boussinesq vessel 

Let us choose the following parameters 



<7l = Ol, (72 = 



" 


—i 


" 




" 





" 


i 








> 7 = 











_ 
















i 



(48) 



Suppose also that the operators C, X, B evolve with respect to t using the same formulas as for the vessel 
with vessel parameters o\, o~2, 7 substituted with ai, ai, 7. More precisely, we define the following vessel: 
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Definition 16. Suppose that the parameters 0-1,02,7 are defined in Definition\7\andai,a 2 ,y are defined 
in (|48p . Then a (regular non-symmetric) Boussinesq vessel 2Js uss,reg is a collection of operators and 
spaces and a rectangle R 

Wbouss,^!) = (C(x,t),A c ,X(x,t),A,B{x,t);a 1 ,a 2 ,y,-/,{x,t),a 1 ,a 2 ,^;ri,C a ;'R), (49) 

where the bounded operators C(x,t) : H —} C 3 , A^,~K(x,t) y A : H — > H, B(x,t) : C 3 — > H and a 3 x 3 
matrix function j*(x,t) satisfy the vessel conditions ([5]), ((6]), ©, @, (|9]) and the following evolutionary 
equations 

d 



B 



m „ = -(ABa 2 + B^)a{\ 
= a^(jC-a 2 CA ( ), 

|-X = Bo 2 C. 

at 



(50) 
(51) 
(52) 



the operator X(x,t) is assumed to be mvertible on the rectangle R. In the case A( = A* and C = B* , we 
call such a vessel symmetric. 

Theorem 17. Suppose that Q3s uss,re 9 is a Boussinesq regular vessel, defined in (|49[) and t(x) — 
det(X(a;o, to) _1 X(:r, t)) is its tau function, defined for arbitrary point (xo,to) £ R- Then the coefficient 
3d 2 

q(x) — — ln(r(a;, t)) satisfies the Boussinesq equation flj on R. 

Proof: Consider the Boussinesq equation (fTJ): 



9 r a 2 i 

From the formula (1371) we find that it can be rewritten as 



-(7rii) xti = Q^2[-2 nxxx ~ 4 (~2^' Kll ' >x ^ 



or integrating with respect to x once and multiplying by — — it is enough to show that: 

(7Tii)« = {■Kii)xxxx + 12(^11)^(7111)^. 



(53) 



The formula (|18|1 . where we substitute the vessel parameters with tilde notation, and substitute x- 
derivative with the t-derivative, implies that 



oi 1 a 2 II n+1 — H n+1 a 2 a 1 1 = (H n )' t —a 1 (7 + a 2 H o- 1 — aiH a 2 )H n + H n 'ya 1 . 
Since 7rii = tr(o2-Ho), we can use ()54[) to obtain: 



(54) 



J^ttii = tx(a 2 H ) = tr(ff 2 J^-ffo) = i(n'i2 ~ "21) = ^ tr ( 



-i 

1 




Ho) 



after cancellations and using (|30[) . So, differentiating again and using (154|1 we obtain: 



dt 2 



d_d_ 

dt dx 



TTii = — — tr( 



-i 
z 




d. 



■tr( 



I 



— i 
z 




= -[7ri27r u J - vrii7r-n - -7Tii7rn + -7rii(-4(7rn) + 7r u ) + -7r n . 



Plugging here (|35p we obtain after cancellations (I53|) . 



□ 
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5.2 Standard construction of a Boussinesq vessel 

We may revise now the standard construction of a vessel, presented in Section [3] in order to show how to 
construct a Boussinesq vessel and hence how to produce a solution of the Boussinesq equation UJ). For 
this we have to continue the process of construction, using the parameters <7i,<72,7 defined in (|48[) : 

1-3 use old steps to construct C(x), X(x), B(x) with the same assumptions, as in Section^ 

4. solve the equation 

^-R(\)B = — (A R(\)B<T2 + R(X)B^/)a^ 1 (55) 
ox 

with the initial condition B(x) resulting in a function B(x,t). We have to require, similarly to the 
construction of B(x) that 

VA £ spec(A) : R(\)B(x,t)a 2 C 3 C D(A), (56) 
VA ^ spec(y4) : R(X)B(x,t)jC 3 C H, (57) 

1. solve for C{x,t) the equation (1511) with initial condition C(x) on D(A(), 

2. solve for X(a;,t) the equation (|52|) with initial X(a;), 

3. define y*(x, t) by (O at all points, where \(x,t) is invertible. 

Meanwhile, we present a weaker, regular case form of a general Theorem [19] on producing solutions of 
©: 

Theorem 18. Suppose that a collection 

(C(x, t),A c ,X(x, t),A, B(x, t); o\, a 2 , 7, 7* (a;, t), a\, a 2 , 7; "H, C 3 ; R) 

is obtained from a regular (i.e. all the initial operators are bounded) 5(A) £ 1Z. by the standard con- 
struction around a point (xo,to). Then this collection is a Boussinesq vessel on a rectangle R, including 

3 d 2 

(xo,to). The coefficient q(x) = -— tt(X (xo,to)X,(x,t)) satisfies the Boussinesq equation (IT]) on R. 

2 <7X 

Proof: Notice that since X(so,t2) is invertible, by a continuity, since all the operators are bounded, it 
will be invertible on a rectangle R, including (ieq, to). Then the Theorem follows from the definitions and 
Theorem [171 □ 

6 Examples of solutions of the Boussinesq equation ([ID 
6.1 Solitons 

If we take H = C, we will obtain solitons. Use the following choice of solutions of the vessel equations 
A=(2i^) 3 , AC=(ifj,) 3 , 



B(x,t) = 



E B -2ictfiE B (2a/u) E B \ , E B = exp(2ia^(x + 2afit)) 



Ec = exp(ia/j,(x — afit)), 



Ec 

C(x,t) — —ia/iEc 
{ptiifEo 
s/3 

2acosh( — fi(x + fit)) E x /./o ./o \ 

x,t) = 2_ , E X = exp ^ M x(2a 2 - 1) + ^fi 2 t(2a 2 + 5) 
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where a — e 3 is the basic third root of 1. Then simple calculations show that 

9^ 2 



q(x) = -- 



2 cosh ( — /j,(x + tfi)) 

is the Boussinesq soliton solution of fl]) (see |SCM73j V 

Another non-trivial soliton, obtained in this form is as follows. 



A = (2i M ) 3 , AC, = (2i M ) 3 , 
B(x,t) 



ki = a2i[i, fi — a 2ifi, 



Eb Ej 



E B 



Ef 



-bE E 



cE f 



b = 



1 



C{x,t) 



X(x,t) = 
and the coefficient q(x,t) 



% 

A 

2ix/j, ^2y/3xfi 



+ e 



64^ 5 



q(x,t) 



, 2 \2 



(58) 



Eb = exp(fcia; — ik-^t) 
E f = exp(/ia; - if it), 



(59) 



6.2 Solutions, belonging to the Schwartz class. 

Suppose that we are given real coefficients q(x),p(x) in the Schwartz class and define iru by formula 
(|37[) . Suppose that the functions Fi,F2,Fz axe solutions of (0) satisfying 



Fi{x,k) =e~ kx (l + 0{l)), 
F 2 {x,k)=e~ akx {l + 0(1)), 
F 3 (x,k)=e- a2kx (l + 0(l)), 



as x — y +oo, 
as x — > +oo, 
as x — > +oo 



2n 2tt 

where k £ = {fc| — < arg(fc) < ~^~}- F° r the coefficients q(x),p(x) € S, where S - the Schwartz 

class of rapidly decreasing functions, the existence of Fi,F2,F 3 is shown in |DTT82) . for example. 

Then a function-matrix ty t (\,x) which is a solution of Q12p can be constructed from solutions 
2, *.,3 of © and is given in view of (J2TJ by (k = v^A) 



f*(\,x) = 



1 

3a 2 



where 



-7T11*.,1 - *'„,! 



*,,2 
-7rn*, l2 - *'»,2 



**,3 
-7T11*.,3 ~ *'*,3 



_ 7T21*„,1 ~ (7Tll*»,l + *'„,!)' 7T 2 1*»,,2 - (tTii*,,2 + *'», 2 )' 7T21*,,3 - (7Tll#«,3 + #',,3)' . 



(60) 



=a 2 (Fi+F 2 + F 3 ), y, >2 = —Fi + T F 2 + T F 3 , #„, 3 = -^Fi + —F 2 + -^F 3 



k 



k 2 ' 



k 2 ' 



k 2 ' 



1G 



As a function of A the matrix-function >3/„(A, x) is analytic in the whole plane except for the cut along the 
negative axis (— oo,0]. Notice that this function is globally bounded as there are asymptotic formulas for 
such solutions appearing in [DTT82| . Actually, the three functions Fi,F2,Fs have the asymptotic be- 
havior exp(a 8 fc:r)(l + 0(fc)), i = 0, 1, 2 as k — > oo for all x £ R and these asymptotics can be differentiated 
infinitely many times. 
So we define 

S(\,x) = **(A,a;)$' 1 (A,a;), 

which becomes 

S(X,x) = <i>4x,X)S(\,0)^ 1 (X,x) 

with &,(x, A) - the fundamental solution of (|12|) equal to identity at x = 0. Since Vt* and $ are analytic 
functions of k in Q3, they will be analytic functions of A, except for the cut along the negative real 
axis, where it will have jumps, but is remained bounded at the infinity. These jumps are expressible by 
parameters Ri, R2, appearing at [DTT82) in a quite general case [DTT821 Theorem 14] of real coefficients 
q(x),p(x). 

Applying the standard construction to S(A, 0) we will obtain a vessel realizing qv(x) and derived from 
it pv(x), which is equal by Lemma [151 to the given q(x): qv(x) = q(x). Evolving further the operators 
with respect to t, we will create a solution of fTJ, coinciding with q(x) at time t — 0. 



6.3 General solutions 



We can obtain more general solutions of {TJ if we choose to work in a more general setting, presented in 
Section Applying the standard vessel construction to a function S(X) £ H, we obtain 

Theorem 19. Suppose that a Boussinesq vessel VGbouss is obtained from 

S(X) =1- (7 a X _1 (A7 - A)- x B oi 
by the standard construction, satisfies the following assumptions (i(x,t) £ M. 2 ): 



B(x,t) 



£ D(A), B(x,i) 



£ D(A), B(x,t) 



£ H. 



(61) 



Then the coefficient q(x, t) is twice differentiable with respect to t and four times differ entiable with respect 
to x on an open subset Q C (R x M)\Z, obtained after removing a closed subset Z C R x R of the points 
(x,i) £ Z, in which X(a;, t) is not invertible. Moreover, the coefficient q(x,t) satisfies the Boussinesq 
equation {TJ on fi. 

Remarks: 1. Notice that in the regular case the operator H(x,t) and hence the tau function t(x, t) — 
det(X _1 (a;o, to)~K(x, t)) are analytic on R x R. So, the zeros of r(x,t) are analytic sets of R x R, but in 
general their structure is complicated and out of the scope of this work. 2. The assumptions Q6ip imply 
the regularity assumptions pSfl. (|43)l . ((56)) . ((57)) . 

Proof: Notice that the regularity assumptions (I42[) . (|43[) and (|56p . (|57|) imply that the operators 
B(x,t) : C 3 -> H, G(x,t) : C 3 -> D(A ( ) are defined for all (x,t) £ R 2 (actually, for C(x) = 
[ ci(x,t) C2(x,t) cs(x,t) ] one has to verify that the three entries ci(x,t), oz{x, t), 03(2;, t) are in 
D(A(;) for all x). The operator 'K(x,t) is the unique solution of ([7]), (|52p with the initial condition Xq. 
Notice that this is an operator of the form Xo + T(x, i) for a trace-class operator T(x, i) and we can 
define r(x,t) = det(X,J x X(a;, t)) (O. Moreover, since X(a;) is a;-differentiable 



T~'(x,t) 
r(x,i) 



tr(X'(x,t)X' 1 (x,t)) = tr(B{x,t)a 2 C{x)X- 1 {x,t)) 
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for all points (x,t) G Z — {(x,i) \ r(x,i) 7^ 0}) = {(x,t) | ~K(x,t) is not invertible}. On the other hand, 

3 d 2 t'(x t) 

In t(t i\ follows fi'nm the HiflFerentiahilitv nf - ' 

2dx 2 

r'(x,t) 



the differentiability of the coefficient q(x) = — — — - \xit{x, t) follows from the differentiability of 



r(x,t) 



So, we investigate the existence of the derivatives for 
Applying the formulas (0), (I16|) we obtain that 



t' (x.t) d r'(x,t) ,t (x.t) K o _i _i / 

Trfe) = Ox^ki + (^fey) = tr(B(x,t)(^ S-T^i V 2 )C(,,t)X W)) 



tr( 



1 

1 




C{x,t)X- 1 {x,t)B(x)), 



which is well-defined on f2. Then, similarly, 

r"'{x,t) _ d r"(x,t) r"(x,t) r'(x,t) 
r(x,t) dx r(x,t) i~(x,t) r(x,t) 



tr( 



1 
0-2 

1 



C(x,t)X- 1 (x,t)B(x)). 



Then, differentiating and plugging ((5J, (I16[) again we can obtain that 

r"'(x,t) 

Since 7Ti2 = [ 1 ] C(x, t)X _1 (a;, t)B(x, t) 

Q. Similarly, 7r2i, 7T23, 7T32 are well defined from (|43|) . Finally, differentiating the last expression, we will 
obtain that 

d 2 r'"(x,t) 



in view of (1431) it is a well-defined function on 



dx 2 r(x,t) 
Here, the derivative of 7Ti2 is as follows 



— 3(-k[ 2 TV21 + 7Tl27r 2 i — 7r 23 — 7r 32 ) 



7Ti 2 = j [ 1 ] C(x,t)X~ 1 (x,t)B(x,t) 



[l ] a 1 j t (x,t)C(x,t)X~ 1 (x,t)B(x,t) 



' " 




1 


) 







" 




1 




-r 








-[ 1 ] C(x,t)X~ 1 (x,t)B(x,t)ja 1 



where we used (0), ()16|) . Thus it is a well-defined function on Q. Similarly for 7r 21 . The derivatives of 
T23, 7T32 are as follows 



71-32 = -[0 1 ] C(x,t)X- 1 (x,t)AB(x,t) 

7T23 = [ 1 ] C{x,t)X~ 1 (x,t)AB(x,t) 
which are well defined by the assumption (|61[) on fl. 



— 7T127T13 + 7T33, 



+ 7ri2(7Tl3 + 7T22 — 7T3l) + 7Tll7r 3 2 + 7T33 , 
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The fact that we can i-differentiate q(x, t) twice follows exactly the same lines, using the parameters 
5i , 02 , 7 and the regularity assumption (I61|l . Using (|54p , the first derivative is 

7Tn = i(ni3 — 7T31 + 7rn(7T21 — ^l'/)), 

which is well-defined on £1. For the second derivatives, it is enough to check that 7T13 — 7T3i, tt2i — 7Ti2 are 
well denned. Using (|54|l again 



iTTia = [ f ] C(x,t)X- 1 (x,t)AB(x,t) 
-m X2 = [10 }C{x,t)XT 1 {x,t)AB(x,t) 



+ 7T217T22 + 7Tll(7T23 + K'm) — 7I~33 

+ 7T127T21 + 7T21 + 7T23 + 7rn(7Tl3 — 7T22 + 7T3l) 



which are well-defined on Q by (|6f p . Similarly for tt2i,tt3i and the proof is finished. 
Remark: Notice that one more derivative of tt' 23 produces a term of the form 



□ 



[ f ]C(x,t)X~ 1 (x,t)AB(x,t) 



which may fail to be a well-defined function, because B(x,t) 



may fail to be at the domain of A. 
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